NASA Contract Report No. 194905 


Reciprocity-Based Experimental 
Determination Of Dynamic Forces 
and Moments: A Feas ibility Study 


Istv£n L. Ver and Michael S. Howe 
BBN Systems and Technologies 


(NASA-CR-194905) RECIPROCITY-BASED 
EXPERIMENTAL DETERMINATION OF 
DYNAMIC FORCES AND MOMENTS: A 
FEASIBILITY STUDY Final Report 
(8BN Systems and Technologies 
Corp.) 72 p 


N94-33023 


Unci as 


G3/71 0009921 




May 1994 









REPORT DOCUMENTATION PAGE 


Form Appro wd 
OW0 Wo. 0704-0188 


Public reporting burden for thb oo lection of Worm# km b Mtimetert to overage 1 hour per reeponee, bduding the tfcne lor reviewing inetructbne, eeercNng emting dele tourcee, 
g^hering end meintelnlng the dele needed, end oorrpieting end reviewing the oofteetbn of Wormetion, Send eorrmente reoerding thie burden eetimete or eny other eepec* of the 
ooHec^on of Wormetion. Including luggeetbne for reducing We burden, to Weehington Headquerter* Servtoee, Directorate tor Inform# Ion (^eratlone end Report*, 1215 Jeffereon Dtvto 
Highway, Sule 1204, AHIngion, VA 22202-4302, and to the Offloe of Mwiagement end Budget P#w work Reduction Project (0704*0183), Weehington, DC 20501 


1. AGENCY USE ONLY ( Lmv* bltnk) 

2. REPORT DATE 

May 1994 

3. REPORT TYPE AND DATES COVERED 

Contractor Report 

4. TITLE AND SUBTITLE 

Reciprocity-Based Experimental Determination of Dynamic Forces and 
Moments: A Feasibility Study 

S. FUNDING NUMBERS 

Contract NAS1-19061 
WU 538-03-14-01 

6. AUTHOH(S) " “ ' ” 

Istvan L. Ver and Michael S. Howe 

7. PERFORMING ORGANIZATION NAME(3) AND AOORESS(ES) 

Georgia Institute of Technology 
Atlanta, GA 30332 

BBN Systems and Technologies (Subcontractor) 
10 Moulton Street 
Cambridge, MA 02138 

«. PERFORMING ORGANIZATION 
REPORT NUMBER 

BBN Report No. 7967 
BBN Project No. 621697 

•. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-0001 

10. SPONSORING / MONrTORINO 
AGENCY REPORT NUMBER 

NASA CR- 194905 


11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Kevin P. Shepherd 

Task 7 Final Report prepared by BBN Systems and Technologies for Georgia Tech under subcontract no. 
A-8612-S1. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Unclassified — Unlimited 
Subject Category — 71 


12b. DISTRIBUTION CODE 


13. ABSTRACT (kUximum 200 word%) 

BBN Systems and Technologies has been tasked by the Georgia Tech Research Center to carry Task 
Assignment No. 7 for the NASA Langley Research Center to explore the feasibility of "In-Sifu Experimental 
Evaluation of the Source Strength of Complex Vibration Sources Utilizing Reciprocity." The task was carried out 
under NASA Contract No. NASI -19061. 


In flight H Is not feasible to connect the vibration sources to their mounting points on the fuselage through force 
gauges to measure dynamic forces and moments directly. However, it is possible to measure the intenor sound 
field or the vibration response caused by these structurebome sound sources at many locations and invoke the 
principle of reciprocity to predict the dynamic forces and moments. The work carried out in the framework of 
Task 7 was directed to explore the feasibility of reciprocity-based measurements of vibration forces and 
moments. 


14. SUBJECT TERMS 

Noise; Vibration; Interior Noise 

IS. NUMBER OF PAGE9 

.Text : .48 .. 

Appendices: 14 

18. PRICE COOE 

17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

It. SECURITY CLASSIFICATION 

20. LIMITATION OF ABSTRACT 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 


Unclassified 

Unclassified 

Unclassified 



NSN 7540-01-280-5500 


Standard Form 298 (Rtv. 2-89) 
Pr*»cft>*d by ANSI Sid. Z39-18 

























NASA Contract Report No. 194905 


TABLE OF CONTENTS 

Page 


REPORT DOCUMENTATION PAGE iii 

TABLE OF CONTENTS ▼ 

LIST OF FIGURES vii 

LIST OF TABLES & 

SUMMARY AND CONCLUSIONS . 1 

1. FORMULATIONS OF RECIPROCITY..,,. 2 

1.1 Reciprocity Prediction of Interior Noise Due to Dynamic 

Force and Moment Excitation 3 

1.2 Experimental Determination of the Excitation Forces and 
Moments Through a Combination of Non-Intrusive 

Direct and Reciprocal Measurements 4 

1.3 Simultaneous Multiple Input 6 

1.4 Advantages and Disadvantages of the Different 

Reciprocity Methods 9 

1 .5 Effect of Redundancy on Prediction Accuracy 11 

2. RECIPROCITY TRANSDUCERS 11 

2.1 Design Goals 11 

2.2 Design of the Reciprocity Transducers 12 

Large Reciprocity Transducer 12 

Checkout 15 

Small Reciprocity Transducer 17 


//, Hi car 

SUNK NOT FILMFC 


PAGE 


- V 


. INTENTIONALLY SUNK 



NASA Contract Report No. 194905 


TABLE OF CONTENTS (cont'd) 

Page 

3. RECIPROCITY EXPERIMENTS 23 

3.1 Test Article 23 

3.1 Two Force Excitation 25 

3.2 Moment Excitation 42 

REFERENCES 48 

APPENDIX A: Sensitivity Analysis of the Redundancy 

on Prediction Accuracy A-l 

APPENDIX B: Construction Drawings for the Large 

Pentatondodecahedron Transducer B-l 

APPENDIX C: Construction Drawings for the 

Small Pentatondodecahedron Transducer C-l 


vi 


NASA Contract Report No. 194905 


LIST OF FIGURES 

Page 

1. Determination of In-Flight Forces F x and F 2 Utilizing Reciprocity 7 

2. Photograph of the Large Pemtatondodecahedron Transducer 13 

3. Wiring Diagram of the 12 Woofers in the Large Pentatondodecahedron 

Transducer 14 

4. Magnitude and Phase of the Electrical Impedance, Z=Ul/L Measured for a 

Large Pentatondodecahedron Transducer 16 

5. Volume Acceleration vs. Excitation Voltage Calibration Curve of a Large 

Pentatondodecahedron Transducer 18 

6. Volume Acceleration vs. Excitation Current Calibration Curve of a Large 

Pentatondodecahedron Transducer 19 

7. Linearity Check; Sound Pressure Level at 1 Meter Minus Loudspeaker 

Excitation Voltage Level Obtained for 1 Volt and 10 Volt Loudspeaker 
Voltage 20 

8. Directionality Check; Transfer Function (Microphone Voltage /Loudspeaker 

Voltage) Measured at 10 Degree Incremental Rotation of the Large 
Transducer 21 

9. Photograph of the Small Pentatondodecahedron Transducer 22 

10. Directivity Check; Transfer Function (Microphone Voltage /Loudspeaker 

Voltage) Measured at 10 Degree Incremental Rotation of the Small 
Transducer 24 

11. Sketch of the Fuselage Model Identifying Excitation Positions 1 and 2 and 

Receiver Positions A, B, C and D 26 

12. Direct Transfer Function vai/Fi and Corresponding Reciprocal Transfer 

Function Uia/Ia 29 


__ NASA Contract Report No. 194905 

FIGURES (cont’d) Page 

13. Direct Transfer Function Ubi/Fi and Corresponding Reciprocal Transfer 

Function vis/iB 30 

14. Direct Transfer Function Uci/Fi and Corresponding Reciprocal Transfer 

Function vie /ic 31 

15. Direct Transfer Function Udi/Fi and Corresponding Reciprocal Transfer 

Function vi o/iD 32 

16. Direct Transfer Function UA 2 /F 2 and Corresponding Reciprocal Transfer 

Function V2 a/1a 33 

17. Direct Transfer Function UB 2 /F 2 and Corresponding Reciprocal Transfer 

Function V 2 b/ib 34 

18. Direct Transfer Function UC 2 /F 2 and Corresponding Reciprocal Transfer 

Function V2c/k: 35 

19. Direct Transfer Function UD 2 /F 2 and Corresponding Reciprocal Transfer 

Function V 2 d/ iD 36 

20. Driving Point Input and Response Characteristics, Excitation Pos. 1 

Transducer Location C 37 

21. Change in Force Spectra due to Successive Reattachments of the 

Shaker/Force Gauge Assembly at Excitation Position 1 38 

22. Comparison of the Directly Measured Transfer Function Uci /Fic with 

Reciprocity Prediction vic/ici 39 

23. Coherence Functions 40 

24. Effect of Averaging on Reciprocity Prediction; Force Excitation at Pos. 1; 

Transducer Locations A, B, C and D 43 

25. Input Accelerance at Excitation Point 2 44 


viii 


NASA Contract Report No. 194905 


FIGURES (cont’d) Page 

26. Change in Force Spectra due to Successive Reattachments of the 

Shaker/Force Gauge Assembly at Excitation Point 2 45 

27. Prediction of Force F 2 Utilizing Reciprocity 46 

28. Comparison of the Direct and Reciprocal Transfer Function Obtained for 

Moment Excitation 47 


LIST OF TABLES 

Table 1: Reciprocity Relationships For Single Point Force And Moment 

Excitation 3 

Table 2: Prediction Of An Unknown Force Of Moment Invoking 

Reciprocity 5 



NASA Contract Report No. 194905 


SUMMARY AND CONCLUSIONS 

BBN Systems and Technologies* (BBN) has been tasked by the Georgia Tech 
Research Center (Georgia Tech) to carry out Task Assignment #7 for the NASA 
Langley Research Center (NASA Langley) to explore the feasibility of "In-Situ 
Experimental Evaluation of the Source Strength of Complex Vibration Sources 
Utilizing Reciprocity. The Task was carried out under NASA Contract No. NAS 
1-19061. 

In flight it is not feasible to connect the vibration sources to their mounting points 
on the fuselage through force gauges to measure dynamic forces and moments 
directly. However, it is possible to measure the interior sound field or the 
vibration response caused by these structureborne sound sources at many 
locations and invoke the principle of reciprocity to predict the dynamic forces 
and moments. The work carried out in the framework of Task 7 was directed to 
explore the feasibility of reciprocity-based measurements of vibration forces and 
moments. The work performed included: 

1 . analytical extension of the principle of reciprocity to deal with 
a moment and multiple forces, 

2. design, manufacture and performance checkout of two 
reciprocity transducers, 

3. assistance to NASA Langley in carrying out specific 
reciprocity experiments utilizing the special transducers, and 

4. analysis of the experimental data obtained by NASA Langley. 

The key results of the work carried out in Task 7 were : 

1 . a pentatondodecahedron-shaped reciprocity transducer that in 
its transmitter mode of operation works as an omni-directional 
loudspeaker and in its receiver mode as an omni-directional 


* a division of Bolt Beranek and Newman Inc. 
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microphone. The transducer is powerful enough to generate 
sufficient signal /noise ratio in measuring the sound-induced 
vibration of a typical fuselage in the reciprocal experiments 
and sensitive enough to yield sufficient signal noise/ ratio in 
sensing the interior sound field in the direct experiments; 

2. a new, formulation of the reciprocity principle (which does not 
require the volume velocity calibration of the transducer) and 
lends itself much better for the in-situ measurements of 
dynamic forces and moments than traditional reciprocity 
methods; 

3. the analysis of limited experimental data collected indicates 
that further improvement in the vibration response 
measurements techniques will be needed to yield sufficient 
accuracy for the reciprocal measurements of moments. The 
needed improvements are within the capabilities of current 
signal processing and transducer technology; 

4. reciprocity techniques for vibroacoustic calibration of aircraft 
fuselages and for non-intrusive, indirect measurement of 
vibration forces and moments are promising to warrant 
further investigations. 

1. FORMULATIONS OF RECIPROCITY 

To predict accurately the structurebome component of the interior noise and to 
optimize measures to reduce it, it is necessary to know the magnitude and phase 
of the dynamic forces, F, and moments, M, acting on the fuselage and to know 
the magnitude and phase of the transfer functions p/F and p/M describing the 
relationship between the sound pressure, p, in the fuselage and the forces and 
moments which cause it. Reciprocity as an experimental tool can be used for 
accomplishing both of these tasks. The ways reciprocity is used are briefly 
discussed below. 
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1.1 Reciprocity Prediction of Interior Noise Due to Dynamic Force and 
Moment Excitation 

Previous work carried out under NASA sponsorship have dealt with the 
reciprocity prediction of aircraft interior noise due to point force inputs [1, 2]. 
Considering moment M=FL as the superposition of two forces, F, of equal 
magnitude and opposing phase acting at a distance L apart, the reciprocity 
relationship for moment excitation also can be derived [3]. These reciprocity 
relationships are compiled in Table 1 below. 

TABLE 1 

RECIPROCITY RELATIONSHIPS FOR SINGLE 
POINT FORCE AND MOMENT EXCITATION 


Force Excitation 

Moment Excitation 


"'-T 



If 


FitF 

M = FL 


p 

h 



L h 



Direct 


Direct 


PT 







1 



L 

Reciprocal 


Reciprocal 
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In Table 1 the symbols are defined as follows: 

F = Point force excitation 

p = Sound pressure response at the receiver location caused by force F or 
Moment M=FL 

a = acceleration response at the former excitation point caused by point sound 
source of volume velocity Q located at the former receiver location; 
acceleration is measured in the same direction as the formerly applied force 

Q = Volume velocity of the point sound source 

j- ** . 

co = 2rrf = angular frequency 

L = distance between the accelerometers measuring the acceleration response 
difference ^ - 5^ (5, - / jcoL is the angular velocity response 

The arrows above the symbols signify that they are vector quantities. 

1.2 Experimental Determination of the Excitation Forces and Moments 
Through a Combination of Non-intrusive Direct and Reciprocal 
Measurements _ 

During the framework of this task we have also developed a new non-intrusive 
measurement method for in-situ determination of a force or moment acting on 
the fuselage structure. Both the new and the originally proposed reciprocity 
relationships are summarized in Table 2. In Table 2 the symbol U oc represents 
the open circuit voltage of the reciprocity transducer operating at its receiver 
mode and i the current through the transducer when operating in its transmitter 
mode. 


PREDICTION OF AN UNKNOWN FORCE 
OR MOMENT INVOKING RECIPROCITY 
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Table 2 contains two prediction methods labeled as "originally proposed" and 
"improved" respectively. The originally proposed method employs an omni- 
directional microphone as a receiver while performing the direct, non-intrusive 
"direct listening" experiment and an omni-directional point sound source of 
known vol um e velocity calibration while performing the "reciprocal sound 
excitation" experiments. As noticed during preliminary experiments, the 
originally proposed method had the disadvantages of requiring (1) volume velocity 
calibration of the sound source, (2) determination of the acoustic center of the 
sound source, (3) exact placement of the acoustic center at the former 
microphone position, (4) omni-directional radiation characteristics of the source. 
These requirements made it unduly difficult to obtain accurate predictions at 
high frequencies. 

To overcome these difficulties we have developed the improved method (see Table 
2 ). Here the same transducer is used for the "direct listening" and the "reciprocal 
sound excitation" measurements; first as a microphone and then as a sound 
source. No acoustic calibration of the transducer is required. The transducer 
does not need to be omni-directional. The open circuit voltage U 0 c and the 
current through the transducer, i, can be measured with high accuracy. 

1.3 Simultaneous Multiple Input 

The analytical basis for measuring simultaneous multiple excitation is illustrated 
in this section by considering the simple case of two dynamic forces Fi and F2 
which act simultaneously on the fuselage structure at engine mounting points 1 
and 2 respectively, as depicted in Fig 1 A. These two unknown vertical forces 
produce the vibration velocity responses V31 and V41 or acoustic response in form 
of open circuit voltage of the reciprocity transducer U3 and U4 at locations 3 and 
4 . These responses are recorded simultaneously in flight. 

Before the flight test or after it, on the ground, the fuselage is excited by forces F3 
and F4 or the acoustic volume within the fuselage is energized by sending 
current i3 and 4 through the reciprocity transducers. The velocity response V13, 
v 2 3, V14 and v 24 is recorded at the former excitation points in the form of the 
transfer functions V13/F3, V23/F3, V14/F4, V24/F4 for secondary force excitation or 
in form of vi3/i3, V23/i3, V14/4, V24/14 as shown in Fig IB and 1 C. 
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Fig. 1. Determination of In-Flight Forces F x and F 1 Utili 2 ing 
Reciprocity 

A. In-flight, record responses v 3 , v 4 or U 3 , U A due to 

unknown forces F x and F 2 . 

B. On ground, apply excitation F 3 or i 3 and record 

responses v 13 and v^. 

C On ground, apply excitation F 4 or i$ and record 
responses v I4 and v^. 
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During these reciprocal experiments the former excitation point 1 and 2 must be 
unrestrained by disconnecting them from the structure from which they receive 
the force excitation during flight (e.g. disconnecting* the engine from its mounts). 
If the reciprocal experiments employ Forces F 3 ,F 4 , as shown on the left side of 
Figs IB and 1C, the two unknown forces F, and F 2 are determined from the 
following equations: 


V3 — Fj (v 3l / Fj ) + Fjfv^ / F 2 ) 
v 4 = (v« / F 1 ) + ^2(^2 / F 2 ) 


(la) 

(lb) 


Since the direct transfer functions in equations la and lb are not known, we 
replace them with their reciprocal counterparts [1,2] given by: 


<V*i)-(Vft) 

(v 32 /F 2 ) = (v 23 /F 3 ) 
(v 41 / F,) = (v 14 /F 4 ) 
(v 42 /F 2 ) = (v 42 /F 4 ) 


(2a) 

(2b) 

(2c) 

(2d) 


This yields: 


v 3 = Fi(v 13 /F 3 ) + F 2 (v 23 /F 3 ) (3a) 

^4 = ^i(v M /F 4 ) + F 2 (v 24 / F 4 ) (3b) 

Equations 3a and 3b can be solved for the sought excitation forces Fi and F2 
yielding: 



(Vj3 / V23) (^14 / ^24) 



F 4 ( v 4 / y 4 ) — F 3 ( v 3 / 1 ) 
(Va»/v M )-(v B /v 13 ) 


(4a) 

(4b) 


Inspection of equations 4a and 4b, shows that both the numerator and 
denominator are in the form of a difference. Consequently, at those particular 

* It is expedient to perform the reciprocal experiments first, before the engine 
is mounted. 
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frequencies where either the numerator or the denominator are very small, the 
possibility of a sizable error may occur, due to measurement inaccuracy. 
Consequently, redundancy (by using a few more shaker excitation locations than 
the number of active excitation locations) is advisable. 


If the reciprocal experiments employ acoustic excitation, as shown on the right 
side of Figs IB and 1C, the force prediction equations takes the form 


p _ *3 (ffi? j ^23 ) ~ *4 (^4 j ^24 ) 

1 [ (vis/vaMYH/yw) 

p _ <4(^4 / ^14 )~ (3(^3 ! ^13) 

2 (V 2 4 / V14 ) — (V^ / V 13 ) 


(5a) 

(5b) 


The analysis presented here was carried out for only two simultaneously acting 
vertical forces. If there are n simultaneously acting forces, each with three 
orthogonal components, x, y, z, then we would need a minimum of n shaker 
locations where the fuselage is excited in the x, y, and z directions (i.e., by 
mounting a small cube-shaped metal block at each location and using this for the 
mounting of the accelerometers and the shakers) or using 3n acoustical excitation 
locations. In both cases, we would have 3n linear equations such as Equations 
3a, 3b or 5a and 5b which can be solved for the 3n or orthogonal components of 
the n simultaneously acting forces. 

1.4 Advantages and Disadvantages of the Different Reciprocity Methods 

As shown in the left and right hand side of Fig. 1, the reciprocal transfer 
functions can be obtained either by force excitation or by acoustical excitation. 
The advantages and disadvantages of each mode of obtaining the reciprocal 
transfer functions are discussed below. 

Force Excitation 

As shown on the left-hand-side of Fig. 1, the reciprocal transfer functions 
V 13 / Fj , V 23 / F 3 / v-j 4 / F 4 and v 24 / F 4 , are obtained by applying the known 
forces F 3 andF 4 by a shaker. These experiments are difficult, cumbersome, and 
frequently not sufficiently accurate. The reasons for this are: 
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a. It takes a long time to mount a shaker and change its position and 
orientation. 

b. Simultaneously attaching many shakers and leaving them inactive (except 
one) changes the dynamics of the structure, and therefore is not 
permissible. 

c. The need to change shaker location or orientation after each set of 
vibration response measurements precludes the automatization of the 
experiment. 

d. Frequently, there is not enough space available in the vicinity of many 
excitation points to accommodate the shaker and force gauge. 

e. Accurate measurement of the force acting on the excitation point is 
difficult because the inertia of the mounting mass (needed for fastening 
the force gauge to the structure) is as large or larger than the input 
impedance of the fuselage. This is usually the case at the frequency of 
structural resonances where input impedance is controlled by structural 
damping. 

Acoustic Excitation 

As shown on the right-hand-side of Fig. I, the reciprocal transfer functions 
v 13 / * 3 / V 23 / is, v 14 / i 4 and v 2 4 / U are obtained by actuating the reciprocity 
transducer as a loudspeaker by forcing the current i 3 or i 4 through it. This 
experiment is simpler than the force excitation experiment because: 

a. no shaker is required 

b. no force measurement is required 

c. the reciprocity transducer can be moved 
more easily than the shaker 
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1.5 Effect of Redundancy on Prediction Accuracy 

A sensitivity analysis was carried out to assess the effect of redundancy* on the 
prediction of unknown excitation forces or moments. The analyses is presented 
in Appendix A. The error in determining excitation strength in the presence of 
noise is represented in the form of normalized mean square error e 2 in 
measuring the true force-velocity products. The analyses indicates that the error 
is proportional to £/(3M)*/2, where M is the number of observation and 
reciprocal excitation locations. M must be greater than the number of unknown 
excitation forces, N. 

2. RECIPROCITY TRANSDUCERS 

This section deals with the design goals, the design, manufacturing and checkout 
of two pentatondodecahedron shaped reciprocity transducers. 

2.1 Design Goals 

The design goals for the reciprocity transducer were: 

a. should work as a transmitter as well as a receiver. 

b. should be linear. 

c. in its transmitting mode of operation should be able to produce 
enough sound power output to produce high sound-induced 
vibration on the fuselage to yield sufficient signal/ noise ratio. 

d. in its receiver mode of operation should have high enough 
sensitivity to yield sufficient signal /noise ratio. 

e. should be small. 

f. should have an omni-directional radiation and receiver 
characteristics. 

g should have a well defined acoustic center. 

* more observation points and more reciprocal excitation locations than the 
number of unknown forces or force components. 
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The actual design, described in the following sections represents a balanced 
compromise between these conflicting design goals. 

2.2 Design of the Reciprocity Transducers 

This section contains a brief description of the design and performance 
checkout of the two reciprocity transducers. 

Large Reciprocity Transducer 

The large pentatondodecahedron-shaped reciprocity transducer is shown in 
Fig. 2. The geometry of the transducer is documented in BBN Drawings 
AA02001A, AA02100A, AA02101A, AA02102A, AA02103A and AA02104A 
presented in Appendix B. 

The 12 moving coil loudspeakers constituting the transducer are 4-inch 
diameter moving coil woofers. Model 206-0408, manufactured by the 
Analog and Digital Systems Corporation (ADS) in Wilmington, 
Massachusetts. These woofers were designed for hi-fi automotive 
applications providing higlTstroke and Tow distortion. Their electro- 
acoustics characteristics are as follows: . ^ T ^ ^ _ 


Coil Resistance 
B1 

Effective Coil Surface Area 
Free Field Resonance Frequency 
Total Moving Mass 
Mech. Suspension Compliance 
Max. Peak Coil Displacement 
Max. Cont. Voltage Allowed 


3 Ohm 

4 Tesla x meter 
5.6 x 10' 3 m 2 
60 Hz 

5 grams 

1300 pm/Newton 
£2.5 mm 
lOVeff 


Figure 3 shows the wiring diagram of the 12 woofers constituting the active part 
of the pentatondodecahedron transducer. The series combination of four 
loudspeaker groups, where each group has 3 loudspeakers wired parallel, yields 
an R=4 Ohm nominal resistance for the transducer. 
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Fig. 2 Photograph of the Pentatondodecahedron Transducer 

OUTSIDE RADIUS: 7 INCH (178mm) 

INSIDE DIAMETER: 5.47 INCH (141.5mm) 




R = 4 ohm 


Figure 3. Wiring Diagram of the 12 Woofers in the 
Pentatondodecahedron Transducer 
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Checkout 

After the transducer housing was manufactured the loudspeakers were installed, 
wired and the completed transducer was checked out for proper polarity, 
electrical impedance, linearity and for the effect of acoustic loading on the 
volume acceleration. 

It was found that all the 12 loudspeakers were moving in phase, the transducer in 
its transmitter mode of operation was found linear up to UL = 5V e ff pure tone 
excitation in the entire frequency range from 100 Hz to 5kHz. Figure 4 shows 
measured frequency dependence of the transducer's electric impedance Z=UL/i' 
As expected, the magnitude of the electric impedance peaks at the 156 Hz 
mechanical resonance frequency of the transducer set by the combined inertia of 
the moving parts of the speakers and stiffness impedance of the air volume of the 
pentatondodecahedron enclosure. At this resonance frequency the voice coil 
movement is vigorous and the back-electric voltage, Blv> induced by the coil 
movement is in such a phase that it tends to reduce current through the coil and 
thereby increase the electric resistance. 

The large pentatondodecahedron monopole transducer is equipped with a 
microphone holder that allows the positioning of a 1/4 in diameter Bruel and 
Kjaer microphone. Model 4135 or 4136 so that the membrane of the microphone 
is at the geometric center of the interior volume. The purpose of this microphone 
was to demonstrate that changes of acoustic loading have negligible effect on the 
volume velocity of the transducer. Preliminary experiments performed at BBN 
showed that change in acoustic loading, caused by the vicinity of hard surfaces, 
had negligible effects on the volume velocity output. Consequently, the interior 
microphone is not needed. We recommend that NASA either should keep the 
1 /4-inch microphone permanently installed to monitor the functionality of the 
transducer or the opening for the microphone should be sealed. 


156 Hz 
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Magnitude and Phase of the Electrical Impedance, 
L /l, Measured for the Pentatondodecahredron 
Transducer 
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Figures 5 and 6 show the volume acceleration calibration of a large transducer of 
nominally identical construction as that delivered to NASA. As shown in Fig. 5 
the transducer, operating in its transmitting mode, produces 1 m /sec to 5 
m 3 /sec 2 volume acceleration* for 1 volt excitation voltage. According to Fig. 6 
the volume acceleration output is about 1 m / sec for one ampere excitation 
current. 

Figure 7 shows the transfer function SPL(1 meter) - dBVL as a function of 
frequency from 50 Hz to 5 kHz, where SPL(1 meter) is the sound pressure level 
(in dB re 20pPa) and dBVL is the excitation voltage level (dB re 1 volt) of the 
loudspeaker. The two curves were obtained for 1 volt and 10 volt loudspeaker 
voltage, respectively. Figure 7 indicates that the large transducer is linear and 
that it produces 90 dB to 100 dB sound pressure level at 1 meter distance for 1 
volt excitation voltage. 


Figure 8 shows the transfer function between microphone voltage and 
loudspeaker excitation voltage obtained by NASA at 1m microphone distance. 
The superimposed transfer functions were obtained by rotating the transducer in 
10 degree increments. The close clustering of the curves in Fig. 7 indicates that 
the transducer has an essentially omni-directional sound radiation characteristics 
up to 1300 Hz. 

Small Reciprocity Transducer 

The small pentatondodecahedron-shaped reciprocity transducer is shown in Fig 
9. The geometry of the transducer is documented in BBN Drawings AA02200, 
AA02201 and AA02202, presented in Appendix C. 

The small transducer has a solid pentatondodecahedron-shaped body with a 
drcumradius of 2.627 inches and an inradius of 2.088 inches. The 12 Type 206- 
0127 tweeters, manufactured by the ADS Corporation, are surface mounted on 
the solid body. The tweeters have a 1 inch diameter, half-cone-shaped 
membrane. Each self contained tweeter, with its magnet and back air volume, is 
mounted on a 2 inch diameter disk which is mounted on each face of the solid 
body. 

* Volume acceleration = j cod where Q.is the volume velocity. 
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Fig. 5. Volume Acceleration vs. Excitation Voltage Calibration Curve of a Large 
Pentatondodecahedron Transducer. 
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Pentatondodecahedron Transducer. 
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Fig. 7. Linearity Check; Sound Pressure Level at 1 Meter Minus Loudspeaker 
Excitation Voltage Level Obtained for 1 Volt and 10 Volt Loudspeaker 
Voltage. 


T.F (dB) 



Fraq. (Hz) 


Fig. 8. Directionality Check; Transfer Function 

(Microphone Voltage/Loudspeaker Voltage) 
Measured at 10 Degree Incremental Rotation of the 
Large Transducer.. . 
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Fig. 9. Photograph of the Small Pentatondodecahedron 
Transducer. 
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The wiring diagram of the small transducer is identical with that of the large 
transducer depicted in Fig. 3. The only difference is that the coil resistance of 
each tweeter is 4 Ohm, yielding a total resistance of 5.33 Ohm. 

Figure 10 shows the transfer function between microphone voltage and the 
excitation voltage of the small transducer obtained by NASA at 1 meter 
microphone distance. The superimposed transfer functions werejobtained by 
rotating the small transducer in 10 degree increments. The data in Fig. 10 
indicate that the small transducer does not have an omni-directional radiation 
pattern. Since the large transducer covers the frequency range of interest from 
100 Hz to 5 kHz, which is wider than anticipated, no attempts were made to 
investigate the reasons for the less than desirable directional characteristics of the 
small transducer. 

3. RECIPROCITY EXPERIMENTS 

Reciprocity experiments with two simultaneously acting forces and moment 
excitation were carried out by the NASA Langley research center on a 5.5 ft 
diameter, 12 ft long stiffened cylinder fuselage model. This section contains a 
brief description of the test article and the test results. 


3.1 Test Article 


The reciprocity experiments were carried out by NASA Langley personnel 
utilizing an existing fuselage model made of carbon fibers embedded in epoxy 
resin. The fuselage model has the following characteristics. 


Length: 

Diameter: 

Skin Thickness: 

Skin Young's Modulus: 

Skin Shear Modulus: 
Skin Poisson Ratio: 
Skin Density: 

Ring Frames: 

Stringers: 

Ends: 


3.6m (12 ft) 

1.68m (5.5 ft) 

1.7xl0- 3 m (0.067 in) 

0.13xl0 12 Pa for longitudinal direction 
O.llxlO 11 Pa for circumferential direction 
0.6xl0 10 Pa 
0.3 

1550 kg/m 3 
10 J-section 

22 hat-section, evenly spaced 
3 layers of 1.25 in. thick particle board 
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Fig. 10. Directivity Check; Transfer Function (Microphone 
Voltage/Loudspeaker Voltage) Measured at 10 
Degree Incremental Rotation of the Small 
Transducer. 
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The mode shapes and resonance frequencies of the test article have been 
previously investigated and the results reported [4]. The first structural 
resonance occurs at 64 Hz and the first acoustic resonance of the interior volume 
at 47 Hz. To provide some acoustic damping a number of foam anechoic wedges 
were placed at the far end of the interior volume. 

3.1 Two Force Excitation 

7" • » * - 

In this experiment two coherent forces were applied to the test fixture by two 
shakers mounted 4 in. apart from each other. For the direct experiments each of 
the forces Fi and F 2 were applied individually, then both forces were applied 
simultaneously. The acoustic response was recorded at four interior positions by 
measuring the open circuit voltage of the large reciprocity transducer placed at 
these four receiver positions A, B, C and D. Excitation and receiver positions are 
identified in Fig 11. 


For the reciprocal experiments the large reciprocity transducer was working in its 
transmitting mode and for each position A, B, C and D the sound-induced 
vibration velocity response was measured by miniature accelerometers placed at 
the former two force excitation positions. 


For the direct and the corresponding reciprocal experiments the following 
transfer functions were recorded. 

Corresponding Reciprocal 

No. Direct Transfer Function Transfer Function 


1 Uai/Fi ViA/iA 

2 Ubi/F 1 viB/iB 

3 Uci/Fl v lC/ic 

4 Udi/Fi vio/iD 

5 UA 2 /F 2 v 2 a/*a 

6 UB 2 /F 2 V2b/1b 

7 UC 2 /F 2 v 2 c/ic 

8 UD 2 /F 2 v 2 d/1d 


9 UA/OF 1 +F 2 ) 

10 UB/(Fi+F 2 ) 

11 UC/(Fi+F 2 ) 

12 UD/(Fi+F2) 


25 


NASA Contract Report No. 194905 



Fig. 11. Sketch of the Fuselage Model Identifying 
Excltatidn Positions 1 and 2 and Receiver 
Positions A, B, C and D. 
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The transfer function pairs 1 through 8 are shown in Figs. 12 through 19. 
Observing Figs. 12 through 15, (which display the direct and reciprocal transfer 
functions between Force Excitation Point 1 and the four Transducer Locations A, 
B, C and D) note that the reciprocal transfer functions vi/i agree reasonably well 
with the directly measured transfer functions U/Fi. However, the transfer 
functions obtained for Force Excitation Point 2, as displayed in Figs. 16 through 
19, yield less satisfactory agreement. 

Figure 20A shows the input accelerance (acceleration/Force) measured at Pos. 1, 
and Fig. 20B the force spectrum measured at the same location. Figure 20C 
depicts the open circuit voltage response of the large transducer at Pos. C in 
response to the above force excitation at Pos. 1. The accelerance peaks at 316 Hz 
indicating a local resonance of the test fixture at this frequency. Below 316 Hz 
the structural response is stiffness controlled and above 316 Hz by the structural 
damping. The force spectrum at Pos. 1 peaks around 180 Hz and has a sharp 
minimum at 316 Hz where the first structural resonance occurs. The open circuit 
voltage peaks around 150 Hz where the transducer has its resonance frequency 
and the excitation force spectrum has a broad peak. The large fluctuations of the 
open circuit voltage response are due to acoustic resonances of the interior 
volume of the fixture and to standing wave patterns. 

Figure 21 illustrates the extent of reproducibility of the excitation force 
measurements. After five successive reattachments of the shaker/ force gauge 
assembly at Excitation Position 1, the transfer function TFu = Fn/Uii between 
the measured force Fjj and the applied shaker voltage Uii was measured. The 
first transfer function TFio = Fio/Uio was chosen as a reference to obtain the four 
normalized transfer functions TFn/TFio, TF 12 /TF 10 , TF 13 /TF 10 and TF 14 /TF 10 . 

Nominally all of these normalized force spectra should be identical. The 
differences seen in Fig. 21 are due to slight changes caused by successive 
reattachment of the shaker/force gauge assembly. The spectral differences in 
Fig. 21 therefore represent the maximum achievable collapse between directly 
measured and reciprocally obtained transfer functions and, correspondingly, it is 
representative for the expected accuracy to predict Force F, from reciprocity 

' .si*-'-- i- • ...... 

measurements. 
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Figure 22 shows the ratio of the directly measured transfer function Uci/Fic and 
the reciprocally predicted transfer function vic/ici- Curve A represents the ratio 
of the magnitudes and Curve B the difference in phase. Differences are the 
largest where either one or both of the coherence functions, shown in Fig. 23, has 
a low value. 
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Fig. 12. Direct Transfer Function vai/Fi and 

Corresponding Reciprocal Transfer Function 
UiA/iA- 
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Direct Transfer Function Udi/Fi and 
Corresponding Reciprocal Transfer Function 
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Fig. 19. Direct Transfer Function Ud 2^2 an d 

Corresponding Reciprocal Transfer Function 
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Fig. 20. Driving Point Input and Response Characteristics, 
Excitation Pos. 1 Transducer Location C. 

A. Input Accelerance at Pos. 1 

B. Force Spectrum at Pos. 1 

C Open circuit voltage Response at Pos. C 
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Fie 21 Change In Exictation Force /Shaker Voltage Transfer 

Function Due To Five Successive Reattachments Of The 
Shaker/Force Guage Assembly At Excitation Position 1 
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Figure 24 shows how averaging can improve the reciprocity prediction of 
excitation force Fi . The top curve shows the ratio of the four individual direct vs. 
reciprocal transfer functions measured for transducer locations A, B, C and D. A 
further frequency averaging (e.g. in 1 /3-octave bands) would result in even 
smaller differences. 

Figure 25 shows the driving point accelerance (acceleration/ force) at Pos. 2. It is 
substantially different from the driving point accelerance at Pos. 1 depicted in 
Figure 20. 

Figure 26 shows the change in force spectra at excitation Pos. 2 due to successive 
reattachments of the shaker/ force gauge assembly. When compared with Fig. 
21, obtained for excitation Pos. 1, note that the repeatability at Pos. 2 is not as 
good as was in Pos. 1. 

Figure 27 shows the results of predicting force F 2 from the reciprocal transfer 
functions. The top graph shows each of the four individual predictions while the 
middle graph shows the average of these four predictions. Just as was the case at 
Pos. 1, (see Fig. 24) the averaging here also improves the prediction accuracy. As 
shown in the bottom graph, a frequency averaging in 1 /3-octave bands yields 
further improvements in prediction accuracy. 

The error in measuring excitation force is larger than the random error 
(exemplified by the less than desirable repeatability) because it also entails 
systematic errors caused by the fact that the impedance of the shaker mounting 
hardware (located between the force sensor and the test object) is not negligibly 
small compared with the point impedance of the test object, especially at 
resonance frequencies of the structure where the input impedance takes very low 
values. Consequently, further refinement of our ability to measure the excitation 
force accurately should be achieved before a reciprocity analyses is attempted to 
predict the magnitude and phase of two coherent excitation forces Fi and F 2 
when they act simultaneously. 


41 


NASA Contract Report No. 194905 


3.2 Moment Excitation 

The moment excitation to the test fixture was applied by attaching a rigid 1.75 in. 
long beam stub perpendicular to its skin and applying a force in a direction 
perpendicular to both the cylinder axis and to the beam stub. This way the 
cantilevered beam stub exerted a combination of a moment and an in-plane force 
on the test fixture's skin. The direct transfer functions were defined as the U/F 0 
where U is the open circuit voltage of the reciprocity transducer and F 0 is the 
force applied to the beam. The corresponding reciprocal transfer function was 
v 0 /i where v 0 is the velocity response at the former force excitation location and i 
is the current driving the reciprocity transducer. Figure 28 shows one of the 
representative pairs of transfer functions. The agreement between the transfer 
function pairs is unsatisfactory. It is believed that the lack of agreement results 
from the fact that in measuring the reciprocal transfer function, the predominant 
motion of the skin of the test fixture is radial and the cross sensitivity of the 
accelerometer is not sufficiently low to enable it to measure the much smaller 
circumferential motion of the beam end caused by the rotational motion of the 
skin. 

It is believed that improvements in measurement techniques and signal 
processing procedures can be found which will enable us to perform reciprocity 
prediction of exiting moments with sufficient accuracy for engineering use. 
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Fig. 27. Prediction of Force F 2 Utilizing Reciprocity 


A. Individual prediction from transducer 

locations A, B, C and D 

B. Average of the four predictions 

C Curve B, frequency averaged in 1/3-octave 
bands 
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§1 Formulation of an analytical scheme 

The Reciprocal Theorem Consider a mechanical system of elastic structures 
immersed in a stationary acoustic medium, as indicated schematically in Figure 
1. Forces , F 2 , . . . , are applied at the points labeled respectively 1, 

2,..., in the figure. These forces produce velocities V a , V 2 , V 3 , . . . , at 
points labeled I, 2, . 


\ F 1 

T F 2 

< 

i 

\ 

N V 2 

tA 


2 

w 

I 

1 

1 


Figure 1 The direct problem 

Consider next a reciprocal problem in which forces * / 2 > / 3 are 

applied at 1, 2, 1, producing velocities v x , v 2 , at points 1, 2, 

(Figure 2). 


s 

' V 1 

i V 2 

At 

\^Z 

i A 

2 

I 

2 

— — t — — 


Figure 2 The reciprocal problem 
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According to the Reciprocal Theorem [1], 


V v > + V v 2 + ••• + A- V > * ••• • 

The theorem remains valid if one or more of the forces are replaced by 
couples, in the following sense. Let a couple Mj be applied at point 3 in the 
direct problem, and let the reciprocal forces / 3 , / 2 , / 3 , . . . , produce an 
angular velocity « 3 at point 3 in the structure, then equation (1) becomes, 

V v i + V V 2 + V W 3 ••• -A* V : 1 + / 2 * V 2 + A- V 3 + • < 2 > 

To take account of all possibilities we shall adopt a notation in which, 
F x , F 2 , Mj , etc, are denoted by the generalized forces P x , P 2 , P 3 , etc, and 
the corresponding velocities and angular velocities v 3 , v 2 , w 3 , etc, are 
replaced by generalized velocities q 3 , q 2 , q 3 , etc. The reciprocal theorem 
may then be expressed in the composite form, 

N M 

I *q 4 - I /i* v i . < 3 > 

i-l i-l 


where N and M are arbitrary. 

Application to source strength evaluation Equation (3) forms the basis for 
the experimental determination of unknown forces P^ arising in a "direct" 
experiment in terms of observed values of the velocities V and of the 
velocities q i when auxiliary, reciprocal experiments are performed with known 
forces / . It is assumed that experimental data are available for the 
velocities V from observations of the direct experiment, and that we are at 
liberty to choose the forces f in the reciprocal experiments in an arbitrary 
manner. In each of these latter experiments all of the f i are zero except one 
component of one of the forces, making a total of 3M experiments in which the 
velocities q A are to be measured. 

To simplify the description of the analysis we introduce the notation: 

The forces P i , i-l, 2, 3 N, are represented by the 3N-dimensional 

vector P i , i-l, 2, 3,...,3N. This means that P 3 - (P l ,P 2 ,P 3 ), P 2 - 
(P A ,P 5 ,P 6 ), etc. 


1 . 



PI 
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2. The velocities , i - 1, 2, 3 M, are represented by the 3M- 

dimensional vector , i - 1, 2, 3.....3M, so that V 1 — (^i>^ 2 ’^ 3 ^’ ^2 " 

<WV» etc - • 

3. In the 3M reciprocal experiments the forces f i are defined by f 1 - 
(/ 1 ,/ 2 ,/ 3 ), / 2 - (/ 4 ./ 5 ,/ 6 ). etc., where in the first reciprocal 
experiment f ^ is the only non-zero component of force, in the second only 
f 2 is non-zero, etc. 

4. In the jth reciprocal experiment the generalized velocities q i i - 1, 2, 

3,...,N, are denoted by the 3N-dimensional vector (i — 1 to 3N) , 

where q x - (q^ ,q J2 ,q J3 ) , q 2 - (q J4 )» etc - 

For each of the 3M reciprocal experiments the reciprocal theorem can be 
applied in the form, 


3N 

l 

i-l 


V P i 




j - 1, 2, 3 3M. 


(4) 


The quantity is a matrix of order 3Mx3N (i.e., 3M rows and 3N columns). 

To calculate the unknown forces P i from this system of equations the number of 
reciprocal experiments 3M must be at least as large as 3N. In practice, 
however, the data relating to the direct and reciprocal velocities V i , q^ 
will be contaminated by noise, and it will be advisable to use an 
overdet enzined system of equations to improve the accuracy of the predictions. 
This means that M > N, and equations (4) must then be solved by the method of 
least squares [2], according to which the solution vector of (4) is chosen 
to minimize the squared residual 


3M f 3K 

ri ■ M-i vP ‘ ■ /jVj ^ 


(5) 


Taking the derivative of this equation with respect to P i and equating 
the result to zero, we obtain the following system of equations for the 
determination of the : 

3M 3N 3M 

l l q 31 q Jt p t - l 1 - l > 2 - 3 3N - (l 

j-1 1-1 j-1 

These equations are called the norms! equations of the system. Observe that 


A-3 
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Z^q is a square matrix of order 3Nx3N, and the system (6) is just 

sufficient to determine the unknown forces Pj , i - 1, 2, 3, 3N, in terms 

of the overspecified 4ata furnished by the 3M reciprocal experiments. In the 
unlikely event that all the data have been obtained without error, the 
solution of (6) will be independent of the value of M when M > N. 


§2 Sensitivity analysis 

It may be assumed that the amplitudes of the reciprocal f orce-velocity 
products /jV are all normalized to unity. To simplify the notation we write 

[3], 


3M 

a ij " l \i ■ *i - A v i - 

J-i 

and define the 3Nx3N matrix to be the inverse of a i . , i.e. p 

3N 


(7) 


£ c ik Q kj - 


k-l 


( 8 ) 


Assume that the respective errors e. , t in the force-velocity products 
y i -■ / i V i , y. - /^V are uncorrelated for i * j, and that the mean square 
errors e 2 , say, are equal, i.e., 


<€ i € " € 2<5 ij * 


(9) 


where the angle brackets < > denote a time average. The corresponding error 
S? x in P i is given by, 


3M 


3P, 


sp ‘ ' A '<■ 


Hence , 


<(SP )!> - l i « « >£ l 1 — 1 - f ‘1 

1 1 ay, 8 5-, j-! U yj J 


3M 


(10) 


( 11 ) 


To calculate 5P i /3y J we use, 


3N 3M 

p i - I c l <? Jk yj ■ 

k-l j-1 


( 12 ) 
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so that, 


3P 


3N 


U i I 

’ k-i c “ " Jk ' 


(13) 


Thus , 

3H 3K 3M 

<(5P i )2>- l l C lk C. £ l < 2 q jk q j£ 
k-1 i-1 j-1 

3N 3N 

■ « 2 Z I c. k c it a kt . 

k-l i-1 

Using equation (8), this finally yields, 


(14) 


<(5P i )2> - < 2 C U . 


(15) 


Example Consider the case of two generalized forces P^, P^ » and write 
equations (4) in the form, 






^3M ^3M J 


V 


V P 2 



y 2 

y 3 

y 3M - 


where A - (q n ,q 21 .... ,q 3M x ) , M - (q 12 .q 22 
(16) can be written, 


(16) 


,q 3M2 ). In shorthand notation, 





The transposed matrix 




(17) 


a 


U 


r A2 A«m ' 
k M 2 . 


(18) 
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It follows that, 


A»y - 

|A|[l-(A.i)2] 

A A A A 

ji-y - (A-;Q(A« y) 

m| (i-(A*i) 2 ] 


- A/| A | , 


P - P/ 


The corresponding error estimates are, 


<(S P 1 ) 2 > 


<(5P 2 ) 2 > 


|A|yti-(A-i)2 


|p|/[l-(A*p) 2 ] 


Provided /[ 1 - (A*/j) 2 ] remains bounded, equations (21) imply that the errors 
decrease like e/(3M) 4 as M becomes large. 

It can frequently happen that the normal equations (6) are ill- 
conditioned. The matrix a., is then very nearly singular (in the two 

A A 

dimensional example above it corresponds to ./( 1- (A*/i) 2 ] being nearly equal to 
2 ero). When this occurs the normal equations can be solved by the method of 
"singular value decomposition", which is described in detail in Reference 3. 
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Construction Drawings for the 
Large Pentatondodecahedron Transducer 
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APPENDIX C 


Construction Drawings for the 
Small Pentatondodecahedron Transducer 
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